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ABSTRACT 

Supernova (SN) rates are potentially powerful diagnostics of metal enrichment and SN physics, 
particularly in galaxy clusters with their deep, metal-retaining potentials and relatively simple star- 
formation histories. We have carried out a survey for supernovae (SNe) in galaxy clusters, at a redshift 
range 0.5 < z < 0.9, using the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. 
We reimaged a sample of 15 clusters that were previously imaged by ACS, thus obtaining two to 
three epochs per cluster, in which we discovered five likely cluster SNe, six possible cluster SNe la, 
two hostless SN candidates, and several background and foreground events. Keck spectra of the host 
galaxies were obtained to establish cluster membership. We conducted detailed efficiency simulations, 
and measured the stellar luminosities of the clusters using Subaru images. We derive a cluster SN rate 
of 0.35 SNus ^°;i2 (statistical) ±0.13 (classification) ±0.01 (systematic) [where SNus = SNe (100 yr 



lOi°Ls,0)"^] and 0.112 SNua/ lo.ogg (statistical) ±0.042 (classification) ±0.005 (systematic) [where 

SNua/ — SNe (100 yr 10^° Mq)~^]. As in previous measurements of cluster SN rates, the uncertainties 
are dominated by small-number statistics. The SN rate in this redshift bin is consistent with the SN 
rate in clusters at lower redshifts (to within the uncertainties), and shows that there is, at most, only 
a slight increase of cluster SN rate with increasing redshift. The low and fairly constant SN la rate 
out to z K, 1 implies that the bulk of the iron mass in clusters was already in place hy z k, 1. The 
recently observed doubling of iron abundances in the intracluster medium between z — \ and 0, if 
real, is likely the result of redistribution of existing iron, rather than new production of iron. 

Subject headings: supernovae: general - galaxies: clusters: general 
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1. INTRODUCTION 

Quantifying the rates and properties of supernovae 
(SNe) in high-redshift galaxy clusters is important for 
several applications. In structure-formation studies, SNe 
play a crucial role in baryonic physics. Their energy de- 
position into the environment is relevant to both galaxy 
formation and star formation. Numerical simulations of 
galaxy formation now include feedback from SN explo- 
sions (e.g., Borgani et al. 2004; Kay et al. 2007; Nagai 
et al. 2007; Scannapieco et al. 2008; see Borgani et al 
2008a, b for reviews), but the efficiency of this feedback 
is unknown. 

In terms of cosmic metal-enrichment history, SNe are 
the sources of iron and other heavy elements that can 
be observed in the intracluster medium (ICM) and are 
detectable through X-ray observations (e.g., Balestra et 
al. 2007; Maughan et al. 2008; de Plaa et al. 2007). 
The abundances of these elements in the ICM depend 
on the integrated history of SN explosions (e.g., Maoz & 
Gal- Yam 2004), as all of the elements produced during 
all stages of cluster formation and evolution must remain 
in the cluster due to its deep potential well. The abun- 
dances also depend on the efficiency with which matter 
is ejected from galaxies into the ICM, whether by SN- 
driven galactic winds (De Young 1978; White 1991; Ren- 
zini 1997; Borgani et al 2008b; Sivanandam et al. 2009), 
by gas stripping due to ram pressure (Gunn & Gott 1972; 
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Mori & Burkert 2000), or by galaxy-galaxy interactions 
(e.g., Clemens et al. 2000). SNe from a diffuse intergalac- 
tic stellar population may also be non-negligible contrib- 
utors to the ICM enrichment (Gal- Yam & Maoz 2000a, 
2000b; Gal- Yam et al. 2003; Lin & Mohr 2004; Torna- 
tore et al. 2007). Measuring the properties and rates of 
SNe of all types in clusters as a function of redshift can 
thus shed light on galaxy and cluster formation. 

Finally, cluster SN rates can provide clues for our un- 
derstanding of Type la supernova (SN la) physics. It is 
widely agreed that SNe la are the thermonuclear explo- 
sions of near-Chandrasekhar-mass carbon-oxygen white 
dwarfs in binary systems. However, the nature of the 
progenitor systems is still not known, and several dif- 
ferent channels have been proposed (see, e.g., Mannucci 
et al. 2008 for a recent overview). One prediction of a 
progenitor scenario that can be tested by observations is 
the delay-time distribution (DTD) between formation of 
a stellar population and the SN la explosion of some of 
its members. In recent years, constraining the DTD has 
been attempted by comparing cosmic star-formation his- 
tory (SFH) to redshift-dependent rates of SNe la in the 
field (e.g.. Gal- Yam & Maoz 2004; Dahlen et al. 2004, 
2008; Cappellaro et al. 2005; Neih et al. 2006; Bot- 
ticella et al. 2008; Poznanski et al. 2007; Kuznetsova 
et al. 2008). A major complication in such measure- 
ments is the observational uncertainty in the SFH (e.g., 
Forster et al. 2006). In a further recent development, 
several studies have found evidence for the coexistence 
of two SN la explosion channels, a "prompt" channel 
that leads to an explosion within ^ 10* yr of the forma- 
tion of a progenitor binary system, and a "delayed" one 
that occurs at least several Gyr after star formation and 
dominates in old stellar environments (Mannucci et al. 
2005, 2006; Scannapieco & Bildsten 2005; Sullivan et al. 
2006b; Totani et al. 2008; Pritchet et al. 2008; Aubourg 
et al. 2008; Raskin et al. 2009; Maoz et al. 2010; Brandt 
et al. 2010; Maoz & Badenes 2010). 

Galaxy clusters form unique environments for DTD 
studies. The stellar population in galaxy clusters is dom- 
inated by old stars in early- type galaxies, particularly 
in the core of the cluster (e.g., Visvanathan & Sandage 
1977; Renzini 2006) with a very small amount of star 
formation taking place, mainly in star-forming galaxies 
at the outskirts of the clusters (e.g., Hansen et al 2009; 
Porter et al. 2008; Bai et al. 2007; Saintonge et al. 2008; 
Loh et al. 2008). The fraction of star- forming galaxies 
increases with redshift (Butcher & Oemler 1978, 1984), 
an effect that is independent of cluster richness (Hansen 
et al 2009). Since clusters have little ongoing star for- 
mation, measuring the redshift-dependent SN la rate in 
clusters can isolate the delayed channel. Moreover, since 
the SFH in clusters is relatively simple compared to field 
galaxies, the DTD that can be deduced from the SN rate 
depends less strongly on the details of the assumed SFH 
(Maoz & Gal- Yam 2004). 

To date, the SN la rate as a function of redshift in 
clusters has not been well measured. Until recently, the 
few existing published rates relied on small numbers of 
detected SNe, and the large uncertainties were domi- 
nated by small-number statistics. An intermediate/high- 
z SN la rate was derived by Gal- Yam et al. (2002) using 
archival Hubble Space Telescope (HST) imaging of 9 clus- 
ters, in which they discovered two or three likely cluster 
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Fig. 1. — Candidate classification. "BG" and "FG" denote back- 
ground and foreground events, respectively. 

SNe. Their measured rates were 0.39l[5;^^ and 0.80tJ]:^g 
SNub aX z — 0.25 and z = 0.9, respectively, where SNus 
denotes SNe (100 yr lO^^Ls,©)"^. These rates corre- 
spond to roughly O.lltJJ:^^ and 0.22t°;^^ SNum, respec- 
tively, where SNum denotes SNe (100 yr IO^^Mq)"! 
(see § [5]) . Based on three cluster SNe and three possi- 
ble cluster SNe, Graham et al. (2008) derived a SN la 
cluster rate at 0.2 < z < 1.0 from the CFHT Supernova 
Legacy Survey (SNLS) of O.lt^;^^ SNua/. Mannucci et 
al. (2008) have reanalyzed the Cappellaro et al. (1999) 
nearby SN sample to derive a local (z < 0.04) SN la clus- 
ter rate, based on 11 SNe, of O.066l^;5520 SNua/, which 
they found to be significantly higher than the correspond- 
ing rate in field elliptical galaxies, 0.019^g'Qgg SNuM. 

Sharon et al. (2007) reported a rate of 0.098tJ^ J^^l^±0.009 
SNuA/ at a shghtly higher redshift, 0.06 < z < 0.19, 
based on the Wise Observatory Optical Transient Sur- 
vey (WOOTS) detection of six cluster SNe la (Gal- Yam 
et al. 2008). The Sloan Digital Sky Survey-II (SDSS-II) 
Supernova Survey has discovered thousands of SN candi- 
dates, and by cross correlation of the confirmed SNe with 
SDSS cluster catalogs (Koester et al. 2007, Miller et al. 
2005) measured cluster rates of 0.060lES{!20-o.ooi SNua/ 
at z = 0.084 and QMStofiit°ofQ2 SNua/ at z = 0.225 
(Dilday et al. 2010). 

Measurements of SN la cluster rates from larger sur- 
veys are ongoing. To name a few, the Palomar Transient 
Factory (Law et al. 2009; Rau et al. 2009) is expected 
to discover thousands of SNe in a footprint of > 8000 
deg^, which will also be used to measure cluster SN rates 
at low redshift. At very high redshifts, the Supernova 
Cosmology Project (SCP, PI Perlmutter) has targeted 
25 clusters with a 219-orbit HST multi-epoch program 
(GO-10496) in which about eight cluster SNe have been 
discovered (Dawson et al. 2009; Melbourne et al. 2007), 
permitting derivation of the cluster SN rate at z > 1. 

Here we present results of an _ffST-based SN survey in 
galaxy clusters at 0.5 < z < 0.9. Throughout the paper 
we assume a flat cosmology, with parameters JIa — 0.7, 
r2,„ = 0.3, and Hq = 70 km s~^ Mpc^^. Magnitudes are 
reported in the Vega-based system unless stated other- 
wise. 

2. HST OBSERVATIONS AND REDUCTIONS 

We were allocated 30 HST orbits to reimage 15 high- 
redshift galaxy clusters, during two observation cycles 
(Programs GO-10493 in Cycle 14 and GO-10793 in Cy- 
cle 15, PI A. Gal- Yam). The clusters were selected to be 
X-ray bright, in the redshift range 0.5 < z < 0.9, and to 
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have been imaged with the Advanced Camera for Surveys 
(ACS) in the past. We also required that archival data 
for each target were non-proprietary at the time of our 
imaging, to ensure prompt detection of SN candidates. 
To enable comparison to archival data, the new images 
were obtained using the same filters as the archival ones, 
either the F814W filter (~ / band) or the F775W fil- 
ter i band). New observations were obtained at the 
same position angles as the archival ones, or rotated by 
90°, 180°, or 270°, to allow maximal overlap between 
the images. The archival data consist of observations 
from several HST programs; see Table [I] which also lists 
the other observations that were used in this survey, and 
the area of overlap between epochs. Exposure times were 
typically one orbit (~ 2000 s) per observing epoch in pro- 
grams GO-10493, GO-10793, GO-9292, GO-9744, and 
GO-9836, and at least two orbits 4500 s) in programs 
GO-9033, GO-9090, GO-9290, GO-9722, and GO-10509. 
The 5ct detection limiting magnitudes for a point source 
are typically ~ 26.4 mag in the / band and ^ 26.9 mag 
in the i band in our single-orbit observations, and ^ 27.5 
mag in the two-orbit archival F-band (F555W) images 
that exist for some of the fields. The precise detection 
probabilities as a function of magnitude for variable point 
sources are measured through simulations, as detailed 
below in § 17.1.11 

We were able to obtain high-quality data for all but 
one cluster in the original list during Cycle 14. The 
missing cluster, MACSJ0025-1222, was withdrawn from 
our survey since it had not yet been successfully ob- 
served (GO-9722) at the time of our survey. It was re- 
placed in our Cycle 15 target list with another cluster, 
SDSSJ1004-I-4112, that had already been observed twice 
by HST (Cycle 12, GO-9744; Cycle 14, GO-10509) and 
fits our selection criteria. During Cycle 15, the halt in 
operations of ACS caused the early termination of our 
program; six of the targets were not imaged. We there- 
fore have three epochs for nine clusters, and only two 
epochs for six clusters. 

Each new epoch was split into four dithered subex- 
posures. The subexposures of each epoch were reduced 
using the standard HST/ ACS pipeline, and combined us- 
ing the Multidrizzle routine (Koekemoer et al. 2002) to 
remove cosmic-ray hits, dead or hot pixels, and other 
artifacts such as trails from satellites that crossed the 
field of view, with a square kernel, pixfrac=1.0 and 
scale=INDEF. Images from consecutive epochs of the 
same field were then aligned and subtracted from each 
other, to form a difference image. Specifically, transient 
candidates were searched for in Epoch I compared to 
Epoch II, in Epoch II compared to Epoch I, and in Epoch 
III compared to Epoch II. We note that point-spread 
function (PSF) matching was not required (see § 17.1.11 
for more details). 

Each difference image was searched by eye promptly 
after the observation, and all transient or variable can- 
didates were noted. Obvious variable stars and known 
active galactic nuclei (AGNs) were removed from the 
follow-up list at this stage. Table [2] lists the remaining 
candidates. Assuming that SN production follows light 
(e.g., Forster & Schawinski 2008), a galaxy was consid- 
ered to be a candidate's host if the candidate was within 
its 2cr isophotal contour, defined as the contour along 
which the galaxy flux per pixel is 2a above the back- 



ground fluctuations, which in practice means the candi- 
date is seen clearly embedded in the galaxy light. From 
the areas enclosed by these isophotal contours, the prob- 
ability for a chance association in a given image is < 2%. 
Three of the 37 candidates do not satisfy this criterion, 
yet they were also chosen as likely candidate-host associ- 
ations due to their small projected distances to their pu- 
tative hosts, less than twice the radius of the 2a isophotal 
contour. In these cases, the probability for a chance as- 
sociation is 6%. When more than one galaxy could be a 
likely host, follow-up spectroscopy was scheduled for the 
additional galaxies as well. 

As detailed below and shown in Figure [TJ among the 
37 candidates, at least 6 are likely cluster events based 
on their host-galaxy redshifts (Figure [2]), 5 of which are 
likely SNe la, and one is likely a core-collapse (CC) SN. 
Two candidates have ambiguous hosts and are possible 
cluster SNe la. Among the other candidates (Figure |3]), 8 
are background transients (BG), 8 are foreground events 
(FG), and 3 proved to be AGNs. Two candidates have 
no apparent host, and are possible cluster SNe la. The 
host galaxies of the remaining 8 candidates are not con- 
flrmed by spectroscopy, and their classification is also 
discussed below. Of these 8 unconfirmed candidates, 4 
are possible cluster SNe la. The remaining 4 were vari- 
able (rather than transient) point sources that were de- 
tected at all epochs and are not clearly associated with 
any host galaxy, and are therefore probably quasars or 
Galactic variable stars. 

3. SUBARU IMAGES 

Our data analysis and derivation of SN rates made use 
of auxiliary ground-based data (see § [71 below) . Eleven 
of the clusters in our survey were imaged as part of a 
study of MACS clusters (Ebeling et al. 2007) using the 
SuprimeCam wide-field imager mounted on the Subaru 
8.2-m telescope at Mauna Kea (Miyazaki et al. 2002), 
covering a field of view of ^ 34' x 27' per cluster. Ta- 
ble [3] lists the clusters for which we have obtained Sub- 
aru data, the observation dates, and the exposure times. 
The 2.5(7 limiting magnitude is typically R = 25.7 mag 
(AB). Details regarding the Subaru data, including data 
reduction and photometry, can be found in Kartaltepe et 
al. (2008). Briefly, the observations took place between 
2000 and 2007 in variable conditions, with seeing rang- 
ing from 0'.'6 to I'.'l. The data were reduced with the 
standard SuprimeCam pipeline (Donovan 2007). Pho- 
tometric zeropoints were derived from overlap with the 
SDSS, or (for clusters outside the SDSS footprint) using 
3 s exposures of nearby SDSS flelds interlaced between 
the cluster observations. The zeropoint uncertainty is 
0.1 mag, and results from the different Subaru and 
SDSS fllter passbands, and from the non-simultaneity of 
the interlaced calibration photometry. 

Object catalogues for each cluster were created us- 
ing SExtractor (Version 2.4.3; Bertin & Arnouts 1996) 
in "dual-image" mode (i.e., detecting objects in one 
image, while performing photometry on another) with 
the i?-band image as the reference detection image. 
Star/galaxy separation was based on the SExtractor pa- 
rameter MU_MAX (peak surface brightness above the back- 
ground level). Since the light distribution of a source 
(e.g., its half-light radius) scales with magnitude, stars 
and other point sources populate a well-defined locus in 
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TABLE 1 

Cluster, Fields 



Cluster 


Z 


Coordinates (J2000) 


tpocn i 


Epoch II (GO-10493) 




Epoch III (GO-10793) 




Filter 






RA 


Dec 


date GO 


date 


area*^ 


At'* 


B 


M 


date area'^ 


At' 


B 


M 




MACSJ2214-1359 


0.503 


22 14 57.34 


-14 00 12.2 


2003/10 9722 


2005/08 


8.34 


100.5 


4.4 


14.2 














F814W 


MACSJ0911+1746 


0.505 


09 11 11.18 


+ 17 46 34.8 


2004/03 9722 


2005/10 


8.44 


93.4 


2.2 


6.9 


2006/12 10.3 


101.2 


2.5 


7.6 


F814W 


MACSJ0257-2325 


0.505 


02 57 08.83 


-23 26 03.3 


2004/01 9722 


2005/08 


8.05 


75.4 


2.1 


6.4 


2006/08 10.5 


99.0 


2.4 


7.5 


F814W 


MS0451. 6-0305" 


0.538 


04 54 10.48 


-03 01 38.5 


2004/01 9836 


2005/07 


7.96 


83.3 


3.0 


9.4 














F814W 


MACSJ1423+2404 


0.543 


14 23 48.60 


+24 04 49.1 


2004/01 9722 


2006/03 


9.13 


94.4 


3.1 


9.6 














F814W 


MACSJ1149+2223 


0.544 


11 49 35.51 


+22 24 04.2 


2004/04 9722 


2006/05 


9.67 


85.0 


5.2 


15.4 














F814W 


MACSJ0717+3745 


0.546 


07 17 32.93 


+37 45 05.4 


2004/04 9722 


2005/10 


8.67 


83.4 


5.7 


16.9 


2006/10 10.5 


72.2 


6.2 


18.5 


F814W 


MS0016.5+1654" 


0.546 


00 18 32.80 


+ 16 26 06.9 


2002/01 9292 


2006/06 


6.65 


56.4 


1.3 


3.7 














F775W 


MACSJ0025-1222 


0.584 


00 25 30.23 


-12 22 43.0 


2004/10 9722 


























F814W 


MACSJ2129-0741 


0.589 


21 29 26.30 


-07 41 26.2 


2003/09 9722 


2005/06 


8.64 


70.1 


1.7 


5.1 














F814W 


MACSJ0647+7015 


0.591 


06 47 49.78 


+70 14 56.4 


2004/12 9722 


2006/02 


8.41 


71.6 


3.9 


10.3 


2006/11 8.54 


74.6 


3.9 


10.4 


F814W 


SDSSJ1004+4112 


0.680 


10 04 34.72 


+41 12 45.0 


2004/04 9744 


2005/12'' 


4.65 


59.0 


2.3 


12.6 


2007/01 9.68 


58.6 


2.3 


12.6 


F814W 


MACSJ0744+3927 


0.697 


07 44 52.58 


+39 27 26.7 


2004/02 9722 


2005/12 


8.43 


58.2 


4.1 


11.1 


2006/12 8.97 


55.6 


4.0 


10.9 


F814W 


MS1054.4-0321 


0.833 


10 57 00.20 


-03 37 27.0 


2002/12 9290 


2006/01 


9.11 


43.3 


3.8 


11.3 


2007/01 10.4 


43.5 


3.9 


11.6 


F775W 


CL0152-1357 


0.835 


01 52 43.00 


-13 57 20.0 


2002/11 9290 


2005/06 


10.0 


32.4 


2.7 


8.0 


2006/09 7.20 


34.8 


2.5 


7.5 


F775W 


CLJ1226. 9+3332 


0.888 


12 26 58.21 


+33 32 49.4 


2003/04 9033 


2006/01 


10.7 


35.8 


1.1 


3.6 


2007/01 10.5 


37.8 


1.7 


5.5 


F814W 



MS0016. 5+1654 and MS0451. 6-0305 arc also MACS clusters, MACSJ0018. 5+1626 and MACS J0454. 1-0300, respectively. 
^ References for cluster redshifts are Ebeling et al. (2007) for the MACS clusters, Ebeling et al. (2001) for CL0152-1357 and CLJ1226. 9+3332, Tran et al. 
(1999) for MS1054. 4-0321, and Oguri et al. (2004) for SDSS1004+41. 
^ At, Lb, and Al of Epoch I are the same as those of Epoch II. 

GO-10509. 

^ Overlapping area between epochs [arcmin^]. The imaging area of ACS is 10.5 arcmin^. 
^ Visibility time [days], 

^ Stellar luminosity within the search area [10"^^ L^^q] (see ^ I7.2t . 
^ Stellar mass within the search area [10^^ M©] (see I7.4t . 



a MU_MAX/MAG_AUTD plane, and can be excluded from the 
catalog (e.g., Bardeau et al. 2005). In addition, objects 
with peaks sharper than the PSF are not real astronom- 
ical objects, and can be flagged as artifacts. These cata- 
logues were used as supplementary information for can- 
didate classification (§[5l), and for measuring the cluster 
stellar luminosities (§ 17. 2p . below. 



4. HOST-GALAXY SPECTROSCOPY 

Nearly all of the SN candidates were too faint to be 
observed spectroscopically, even at the time of discov- 
ery, and even more so when pre-allocated ground-based 
follow-up observing time arrived. Cluster membership 
was therefore established through spectroscopy of the 
host galaxy. In Table [5J we summarize the spectro- 
scopic information acquired for each of the candidates. 
Follow-up spectra were obtained primarily using the 10- 
m Keck telescopes in Hawaii, either with the Low Reso- 
lution Imaging Spectrometer (LRIS; Oke et al. 1995) or 
the Deep Imaging Multi-Object Spectrograph (DEIMOS; 
Faber et al. 2003), in longslit or multislit configuration. 
The target lists were selected from the HST image, and 
supplemented by other interesting objects in the field, 
such as gravitationally lensed galaxies and cluster mem- 
bers. For DEIMOS multislit spectroscopy, we filled the 
field of view with additional targets drawn from wide- 
field imaging with Subaru or SDSS. Table [5] lists the 
dates, instrument, and observers of each spectroscopic 
observation. 

Multislit observations were reduced and analyzed as 
follows. After standard bias and fiatfield corrections, 
we combined all the observations of the same field into 
a deep, cosmic-ray-cleaned two-dimensional spectrum. 
We compared each multislit spectrum with known night- 
sky lines, and calibrated the wavelength range using the 



IRAF" tasks IDENTIFY, FITCOORDS, and TRANSFORM. Fi- 
nally, we used the IRAF task APALL to subtract the 
background, trace the continuum of the object, and ex- 
tract the one-dimensional spectrum from the calibrated 
image^^. The resulting spectrum was rebinned and 
searched for common galaxy emission and absorption 
lines. Longslit data were reduced in a similar manner, 
and in addition were also flux calibrated using spectra 
of standard stars from the same observing nights. This 
process allows, in principle, a comparison of the spectral 
shape of the object with those of template spectra, and 
derivation of a redshift even in the absence of emission 
or absorption features. In practice, the continuum signal 
from the host galaxies was generally weak, and the red- 
shift determination is based on emission or absorption 
lines. 

5. CANDIDATE CLASSIFICATION 

As our transient candidates are not spectroscopically 
confirmed, we must consider the possibility that some 
or all of them are not SNe. The survey discussed in 
this paper is similar to the one reported by Gal- Yam et 
al. (2002), which was based on archival HST images of 
galaxy clusters, obtained with the Wide Field Planetary 
Camera 2 (WFPC2). While our survey is superior in res- 
olution to the WFPC2 survey, and has a more uniform 
observation scheme, both surveys reach a similar depth, 
the search methods and efficiency simulations are done in 
the same manner, and the time span between epochs is 
similar. As discussed by Gal- Yam et al. (2002), transient 

IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, which are op- 
erated by AURA, Inc., under cooperative agreement with the Na- 
tional Science Foundation. 

Some of the spectra were reduced using tools developed in the 
MATLAB environment (Ofek et al. 2006). 



TABLE 2 
SN candidates" 





Source 






Discovery | 




Host 1 


Classification 


Name 


RA (h m s) 


Dec (° ' ") 


mag^ 


epoch R A (h m s) 


Dec (° ' ") 


mag Morphology Spec z (R,ef)^ 





MACSJ0257-2325-candl 02:57:08.479 -23:24:24.38 25.4 Jan 2004 

MACSJ0257-2325-cand2 02:57:12.425 -23:27:05.08 26.8 Jan 2004 

MACSJ0257-2325-cand3 02:57:10.005 -23:27:14.52 25.8 Aug 2005 

MACSJ0257-2325-cand4 02:57:11.476 -23:27:19.72 26.5 Jan 2004 

MACSJ0257-2325-cand5 02:57:06.904 -23:27:46.79 26.3 Jan 2004 



02:57: 
02:57: 
02:57: 
02:57: 
02:57: 



08.507 
12.425 
09.983 

11.464 
06.947 



-23:24:24.55 
-23:27:05.08 
-23:27:13.54 
-23:27:19.45 
-23:27:46.44 



22.9 
>27 
22.3 
23.6 
20.9 



Irr 

unresolved 
spiral 

E 

spiral 



spiral/Irr 
spiral 
spiral 
unclear 



0.3294 (07) 
N/A 

0.73 (FFl) 

0.4039 (07) 



0.619 (08) 
0.495 (08) 
0.365 (08) 



Foreground 

Possible cluster SN la 

Background 

Possible cluster SN la 

Foreground 



pi 



p 

A 

A 
o 

CO 

o 

O 



MACSJ0647+7015-candl 06:47:30.004 +70:14:54.53 25.8 Feb 2006 

MACSJ0647+7015-cand2 06:47:38.030 +70:16:16.72 26.8 Dec 2004 

MACSJ0647+7015-cand3 06:47:49.824 +70:15:31.19 24.6 Nov 2006 

MACSJ0647+7015-cand4 06:47:59.661 +70:15:19.51 25.2 Nov 2006 



06:47: 
06:47: 
06:47: 
06:47: 



29.79 
38.024 
49.901 
59.626 



+70:14: 
+70:16: 
+70:15: 
+70:15: 



54.12 
16.05 
30.40 
19.37 



22.5 
22.2 
18.6 
23.7 



Background 
Foreground 
Foreground 
Possible cluster SN la 



0.55 (FF2) 
0.4915 (MS) 
0.538 (MS) 
2.084 (FSP) 



MACSJ0717+3745-candl 07:17:38.90 +37:45:20.8 26.3 Apr 2004 

MACSJ0717+3745-cand2 07:17:31.444 +37:44:36.12 24.4 Oct 2006 

MACSJ0717+3745-cand3 07:17:41.444 +37:44:10.54 24.3 Oct 2006 

MACSJ0717+3745-cand4 07:17:40.548 +37:45:06.60 23.3 Oct 2006 

MACSJ0911+1746-cand2 09:11:16.407 +17:47:40.04 24.9 Oct 2005 

MACSJ0911 + 1746-cand4 09:11:16.959 +17:46:48.42 25.1 Dec 2006 

MACSJ0911 + 1746-cand5 09:11:19.290 +17:46:09.27 25.2 Dec 2006 



07:17: 
07:17: 
07:17: 
07:17: 



38.855 
31.519 
41.444 
40.548 



+37:45: 
+37:44: 
+37:44: 
+37:45: 



20.07 
37.58 
10.54 
06.60 



21.7 
20.0 
20.6 
21.6 



spiral 

E 

E 

compact 



0.88 (E) 

0.842 (FF4) 



Cluster, CC 
Foreground 
Cluster SN la 
Background AGN 



09:11:16.407 +17:47:40.04 21.6 



09:11:19.269 +17:46:10.07 21.0 



SO 

No host 
Spiral 



0.553 (FF5) 



Background AGN 

Hostless, possible cluster SN la 

Background 



MACS Jl 149+2223-candl 
MACS Jl 149+2223-cand2 



11:49: 
11:49: 



33.20 
35.734 



+22:24:29.90 
+22:22:18.34 



23.9 
24.5 



Apr 2004 
May 2006 



11:49:33.147 
11:49:35.726 



+22:24:30.35 
+22:22:18.41 



20.0 
23.6 



E 
E 



Cluster SN la 
Possible cluster SN la 



0.87 (GY) 



MACSJ2129-0741-candl 
MACSJ2129-0741-cand2 



21:29: 
21:29: 



28.356 
24.993 



-07:41:34.44 
-07:42:22.76 



25.0 
24.4 



Sep 2003 
Sep 2003 



21:29:28.321 -7:41:34.77 



23.1 



Irr 

No host 



0.582 (FFl) 
0.503 (GYS) 



Background 

Hostless, possible cluster SN la 



MACSJ2214-1359-candl 
MACS J22 14- 1359-cand2 



22:14: 
22:14: 



59.353 
57.012 



-13:58:15.28 
-14:00:11.49 



26.0 
23.7 



Aug 2005 
Oct 2003 



22:14: 
22:14: 



59.382 
57.012 



-13:58:14.92 
-14:00:11.49 



20.9 
19.6 



Spiral 
E 



Background 
Cluster SN la 



0.59 (GYS) 
0.9009 (E) 



CLJ1226.9+3332-candl 
CLJ1226.9+3332-cand2 



12:26: 
12:26: 



54.552 
55.503 



+33:33:56.38 
+33:32:12.42 



25.8 
23.7 



24.9 (i) 
23.5 (i) 



Jan 2006 
Apr 2003 



12:26: 
12:26: 



54.508 

55.588 



+33:33:56.17 
+33:32:12.78 



20.5 
20.2 



Spiral 
SO 



Foreground 

Possible cluster SN la 



19.8 (i) 
21.3 (i) 



0.230 (FF5) 
0.8335 (VD) 



MS1054.4-0321-candl 

MS1054.4-0321-cand2 



10:56: 
10:56: 



56.262 
57.862 



-03:37:51.04 
-03:37:47.77 



Dec 2002 
Dec 2002 



10:56: 
10:56: 



56.287 
57.862 



-03:37:51.59 
-03:37:47.77 



Spiral 
SO 



Foreground 
Cluster SN la 



MS0451.6-0305-candl 



04:54:09.968 -03:00:28.18 



25.5 



Jul 2005 



04:54:10.054 -03:00:27.76 



18.3 



CL0152-1357-candl 01:52:43.099 -13:55:19.76 24.3 (i) Jun 2005 

CL0152-1357-cand2 01:52:37.99 -13:56:25.69 25.2 (i) Jun 2005 

CL0152-1357-cand3 01:52:46.217 -13:58:03.93 24.1 (i) Sep 2006 



23.2 (i) 
21.5 (i) 
23.9 (i) 



Spiral 



0.16 (FF3a) 



Foreground 



1.27 (GY) 
1.12 (FS) 
0.839 (FS) 



01:52: 
01:52: 
01:52: 



43.099 

38.06 

46.217 



-13:55:19.76 

-13:56:25.5 

-13:58:03.93 



compact 

Spiral 

unclear 



Background AGN 
Background 
Cluster SN la 



SDSSJ1004+4112-candl 10:04:33.086 +41:12:31.20 25.2 Dec 2005 

SDSSJ1004+4112-cand2 10:04:30.601 +41:14:10.66 25.9 Dec 2004 

SDSSJ1004+4112-cand3 10:04:31.007 +41:14:13.59 26.3 Jan 2007 



10:04: 
10:04: 
10:04: 



33.075 
30.645 
31.047 



+41:12:30.34 
+41:14:10.67 
+41:14:13.10 



22.3 



21.6 



spiral 



0.753 (FF3b) 
■ • ■ (FF3b) 
0.75 (07) 



Background 

Possible cluster SN la 

Background 



°- Coordinates are J2000. 

^ Magnitudes are observed / band, unless indicated. 

" Redshifts are based on the following observations: (FFl) Keck/LRIS 1" longslit, Filippenko and Foley, 2005 Dec 3; (FF2) Keck/DEIMOS 0.9" longslit, Filippenko and Foley, 2005 Dec 1; 
(FF3) Keck/DEIMOS (a) 0.9" longslit, (b) multislit, Filippenko and Foley, 2005 Dec 31; (FF4) Keck/LRIS l" longslit, Filippenko and Foley, 2006 Dec 20; (FF5) Keck/DEIMOS multislit, 
Filippenko and Foley, 2007 Feb 16; (FF6) Keck/LRIS multislit, Filippenko and Foley, 2007 Jan 12; (FS) Keck/LRIS l" longslit, Filippenko and Silverman, 2007 Nov 12; (FSP) Keck/LRIS 
1" longslit, Filippenko, Silverman, and Poznanski, 2008 Apr 26. (GY) Keck/LRIS, l" longslit. Gal- Yam, 2005 Aug 1; (VD) Keck/LRIS, 1"2 longslit. Van Dokkum at al. (2000); (E) 
Keck/DEIMOS, multislit, Ebeling, 2006 Jan; (GYS) Keck/LRIS multislit. Gal- Yam and Sharon, 2007 Jul 16-17; (MS) Keck/LRIS, l" longslit, M. Sullivan, 2006 Nov 22; (07) Keck/LRIS 
multislit, Ofek, 2007 Jan 22; (08) Keck/LRIS multislit, Ofek, 2008 Jan 4. 
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TABLE 3 
Subaru Imaging Data 



Cluster Name" 




UT Date Observed [mm/yy] 




Exposure time [si 






B 


V 


R 


/ 


z' 


B 


V 


R 


/ 


z' 


MACSJ0257 


09/05 


12/02 


12/00 


12/07 


12/02 


1440 


2160 


5280 


2400 


2700 


MACSJ0454'' 


11/05 


11/05 


03/05 


12/01 


12/06 


1440 


2160 


3240 


2160 


1440 


MACSJ0647 


11/05 


09/03 


02/04 


02/04 


02/04 


1440 


2160 


2880 


2880 


2160 


MACSJ0717 


02/04 


12/02 


12/00 


12/00 


12/02 


1440 


2160 


2880 


1440 


1620 


MACSJ0744 


02/04 


04/03 


12/02 


04/03 


04/03 


1440 


1440 


2880 


3240 


2160 


MACSJ0911 


11/05 


04/03 


04/03 


04/03 


04/03 


2880 


2160 


2880 


1200 


1620 


MACSJ1149 


12/06 


04/03 


04/03 


12/00 


04/03 


1440 


2160 


2880 


1200 


1620 


MACSJ1423 


07/03 


06/02 


06/02 


06/02 


06/02 


1920 


2160 


2400 


2160 


1440 


MACSJ2129 


07/03 


06/02 


06/01 


06/01 


06/02 


2880 


1440 


2880 


2880 


1440 


MACSJ2214 


11/05 


09/03 


07/03 


07/04 


07/04 


1440 


2160 


2880 


2160 


1620 


CL1226 


05/06 


06/02 


12/00 


12/00 


04/03 


2160 


2160 


2880 


1920 


1080 



^ Full cluster names are listed in Table [Tl 

MACSJ0454. 1-0300 is also named MS0451. 6-0305. 



TABLE 4 

Cluster-member SN candidates 



Name 




M% 


R" 


Classification 


Note 


MACSJ0257-2325-cand2 


-14 


7 




0.48 


possible 


Unknown z 


MACSJ0257-2325-cand4 


-15 







0.51 


possible 


Unknown z, early red-sequence hosf^ 


MACSJ0647-l-7015-cand4 


-16 


6 




0.37 


possible 


Unknown z 


MACSJ0717+3745-candl 


-15 


3 




0.56 


not SN la 


Cluster z, likely CC SN 


MACSJ0717+3745-cand3 


-17 


4 




0.69 


likely 


Cluster z, early host 


MACSJ0911+1746-cand4 


-16 


4 




0.49 


possible 


Hostless 


MACSJ1149-|-2223-candl 


-17 


8 




0.32 


likely 


Cluster z, early host 


MACSJ1149-|-2223-cand2 


-17 


2 




0.60 


possible 


Unknown z, early host 


MACSJ2129-0741-cand2 


-17 


5 




0.38 


possible 


Hostless 


MACSJ2214-1359-cand2 


-17 


8 




0.02 


likely 


Cluster z, early host 


CLJ1226.9+3332-cand2 






-18.6 


0.38 


likely 


Cluster z, early host 


MS1054.4-0321-cand2 






-18.9 


0.26 


likely 


Cluster z, early host 


CL0152-1357-cand3 






-18.4 


0.44 


likely 


Cluster z, brightness suggests SN la 


SDSS J1004+41 12-cand2 


-15 


8 




0.67 


possible 


One of the possible hosts is at cluster z 



^ Absolute magnitude at the time of detection, K-corrected to rest-frame V for candidates at ^ < 0.7 and to rest- 
frame B for candidates at ^ > 0.7. In cases of unknown redshift, the SN candidate is assumed to be at the cluster 
redshift. 

^ Projected distance from brightest cluster galaxy, in Mpc. 

^ The host has the same i ~ z' color as the cluster red-sequence galaxies; see Figure 
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Fig. 2. — Discovery images of the possible and likely cluster- 
member SN candidates. In each row of three thumbnails, we show 
the first and second epochs in the left and middle frames, respec- 
tively. In the right panel we show the difference, as a subtraction 
of the earlier epoch from the later epoch. The frames are each 4'.'4 
wide, and centered on the transient candidate. 

events that could, in principle, mimic SNe include solar- 
system objects, variable stars in our Galaxy or in other 
galaxies, AGNs, or gamma-ray burst (GRB) afterglows. 
As Gal- Yam et al. (2002) argue, most of these tran- 
sients cannot be confused with SNe. First, all but two of 
the candidates are clearly associated with galaxies. The 



Fig. 2. — Continued. 

probability for a chance association is small, < 2 — 6% 
(see § [2]) . The proper motions of asteroids or Kuiper-belt 
objects would have been detected in our long-exposure 
images. With the exposure span within one HST or- 
bit (typically ^ 45 min), we can detect proper motions 
greater than ^ 0'.'015 hr~^. At the dates and coordinates 
of our observations, the parallax of a Kuiper-belt object 
at 50 AU due to the Earth's motion is more than 0'.'4, 
well above our detection limit. Variable stars in other 
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Fig. 3. — Same as Figure|2] for non-cluster candidates. 

galaxies are too faint to mimic SNe. 

Variable stars in our own Galaxy, to be undetected in 
one of the epochs, would have to be distant. For example, 
an M5 flare star with absolute magnitude M/ = 9 mag 
in its quiescent state (Allen 1973) would have to be at a 
distance > 40 kpc to be fainter than our I = 27 mag 
detection limit. The local density of halo stars with 
0.09 < M/Mq < 0.71 is of order lO"* pc'^ (Gould et al. 
1998). A number of recent studies (see, e.g., Cignoni et 




MACS2214ca™i1 



4 



ilOMcsndl 

I 4 




Fig. 3. — Continued. 



al. 2007, and references therein) find an approximately 
r^^ profile for the outer halo, and thus the stellar den- 
sity at 40 kpc is 10~4(8.5 kpc/40 kpc)^ w 10"*^ pc'^. 
Integrating this density from 40 kpc to, say, 80 kpc gives 
a surface density of ~ 10^ deg~^. There could thus be 
several hundred outer-halo late-type stars in the 15 ACS 
fields. Based on this, the two hostless candidates we 
have found could, in principle, be distant Galactic M 
stars that had flared into visibility on one epoch. How- 
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Fig. 3. — Continued. 



ever, the flare would have had to be by > 2 — 2.5 mag to 
bring the object from below the detection limit, / = 27 
mag, to the observed brightnesses of the two hostless 
transients, / = 25.1 and / — 24.4 mag, respectively. Al- 
though stellar flares can be bright in the ultraviolet, and 
reach amplitudes of AU ~ 5 mag, their amplitudes in 
red optical bands are much smaller (e.g., Eason et al. 
1992; AUred et al. 2006; Zhilyaev et a. 2007), gener- 
ally AI < 1 mag. Furthermore, a recent SDSS-based 
study by Kowalski et al. (2009) of M-dwarf flaring fre- 
quency and magnitude finds that the flaring duty cycle 
is strongly correlated with Galactic height. Beyond 300 
pc above the disk, only 10~^ of their individual observa- 
tions, which have cadences of several days, catch an M 
star in its flaring phase. Thus, apart from the implau- 
sibility of a large- amplitude /-band flare, there is only a 
~ 10~^ probability to start with, that any of the ~ 100 
outer-halo M stars in the ACS flelds would be caught 
flaring in our observations. It is therefore highly unlikely 
that the hostless candidates are optical flares of Galactic 
stars. 

For candidates that are in the centers of their asso- 
ciated hosts, an AGN nature would have been revealed 
by our spectroscopy. Since GRBs are often associated 



with some CC SNe (Galama et al. 1998; Hjorth et al. 
2003; Stanek et al. 2003; Malesani et al. 2004; Plan et 
al. 2006; see Woosley & Bloom 2006 for a review), we 
argue that in addition to being unlikely, they would have 
affected only our classiflcation as SNe la or CC SNe, not 
the identiflcation as SNe. We thus conclude that essen- 
tially all of our transient candidates are bona flde SNe. 

6. CLASSIFICATION OF CLUSTER SN CANDIDATES 

Cluster membership of a SN candidate was decided 
according to the SN host-galaxy redshift and the cluster 
velocity dispersion (Ebeling et al. 2007). We classified 
SN candidate hosts as cluster members if their redshift 
indicated that their velocity is within 2a of the cluster re- 
cession velocity, where a is the velocity dispersion of the 
cluster. We note that the velocity dispersion of galaxies 
in clusters as massive as those in our sample is relatively 
high; the mean value for MACS clusters in our sample is 
~ 1300 km s"i (Ebeling et al. 2007). 

Next, although we have no spectroscopic classification, 
some clues about the type of the SN candidates exist. 
Since CC SNe are exceedingly rare in non-star-forming 
environments (e.g., Hakobyan et al. 2008), we classify 
all non-AGN candidates in cluster early-type galaxies 
as SNe la. The high resolution of ACS allows classi- 
fication of resolved host galaxies through morphology, 
which can be supported by color information where 
available. In cases where the morphological type is 
unclear, we compare the location of the host in color- 
magnitude space with that of the cluster red sequence 
(Figures H] and E]). For the MACS clusters, we use 
photometric catalogs of multiband Subaru images (see 
§ El for details). Multiband archival HST/ ACS data 
are also available for some of the clusters. The clas- 
sification of each individual candidate is described below. 



6.1. Candidates with Cluster- Member Hosts 

Likely cluster SNe la. — Based on ACS morphology, 
colors, and Subaru photometry, four of the cluster 
candidate hosts are undoubtedly early-type galaxies: 
MACS0717 cand3, MACS1149 candl, MACS2214 cand2, 
and MS1054 cand2. In the three cases where the candi- 
date is not clearly separated from the galaxy core, we 
confirmed via spectroscopy that it is not an AGN. We 
therefore classify these events as likely cluster SNe la. 

CL0152 cand3, likely cluster SN la. — This candidate is 
detected in the second epoch with / « 24.1 mag, slightly 
offset from the center of a cluster galaxy. Although 
the galaxy is bluer than the cluster red sequence, at 
z = 0.8391 the SN absolute magnitude, Mb « -18.4, is 
strongly suggestive of a SN la, given the scarcity of over- 
luminous CC SNe (less than 10% of CC SNe at maximum 
brightness are as bright in a volume-limited sample; Li 
et al. 2010, in preparation). 

MACS0717 candl. — The SN candidate appears in the 
archival epoch, in a spiral galaxy at the cluster redshift, 
and was observed in two filters, V and /. The observed 
V — I color {r^ 1.6 mag) and magnitude (/ « 26.3, corre- 
sponding to My « —15.3 mag at z = 0.548) are consis- 
tent with a CC SN. Young SNe la would be too bright 
(expected to be < 23 mag), and older SNe la would be 
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m0257 z=0.506 



m0647 z=0.58'l 



m2214z=0.504 



m0717z=0.548 



m0454 z=0.55 



m0911 z=0.504 



5 , 
O \ 



m2129 z=0.57 



Fig. 4. — Color-magnitude diagrams for the galaxies in the cluster 
fields with SN candidates for which we have Subaru data. In each 
plot, the brightest cluster galaxy is marked with a black square. A 
dashed line indicates the color of an elliptical galaxy at the cluster 
redshift. Host galaxies are marked with circles (see legend). 



CD152 z=0.831 



Ml 054 z=0.8 
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Fig. 5. — Same as Figure 4, but for cluster fields for which we 
do not have Subaru data. The magnitudes are measured in the 
available ACS bands in the archive. 

too red. We therefore argue that this candidate is not a 
SN la. 

6.2. Candidates with Ambiguous Hosts 

CL1226 cand2 appears in the archival epoch, 1'.'12 from 
the core of an SO galaxy at the cluster redshift. The 
V — I color (~ 2.1 mag) and magnitude (/ ~ 23.7 mag) 
are consistent with a SN la, 9-10 days after maximum 
brightness. However, the candidate is also positioned 
close to a blue emission knot that is possibly a part of 



a disrupted spiral galaxy at unknown redshift, centered 
~ 5" NNW of the cluster member. For SDSS1004 cand2, 
the results of the spectroscopy are ambiguous, with indi- 
cations that the SN may be in a cluster member galaxy. 
However, it is also on the outskirts of a foreground spiral 
galaxy at z = 0.27. Adopting a conservative approach, 
we therefore classify these candidates as "Possible" clus- 
ter SNe la. 

6.3. Hostless Candidates 

Two of the SN candidates, MACS0911 cand4 and 
MACS2129 cand2, have no detectable host at the lim- 
iting magnitude of the coadded ACS images. If these 
are indeed cluster events, they occurred 30 kpc and 
20 kpc away from any galaxy, respectively. The abso- 
lute magnitudes of these candidates. My = —16.4 and 
My — —17.5, are consistent with SNe la at these clus- 
ters redshifts. Their projected distances from the bright- 
est cluster galaxies (BCGs) are 0.49 Mpc and 0.38 Mpc, 
respectively, at the clusters redshifts. The V — I color 
of MACS0911 cand4 does not rule out any SN type at 
the cluster redshift. These, if intergalactic cluster events 
(Gal- Yam et al. 2002), are most probably SNe la, since 
intergalactic CC SNe should be very rare; there is little or 
no star formation in the intergalactic medium (Gal- Yam 
et al. 2003), and even if a CC SN progenitor is ejected 
from its host galaxy, it cannot move far before exploding 
since it is short lived. The progenitors of SNe la, on the 
other hand, have ample time to reach 30-70 kpc from 
their host galaxy before explosion. 

6.4. Candidates Without Measured Host Redshift 

We do not have secure redshift measurements for four 
of the SN candidate hosts. Spectroscopy was attempted 
(unless indicated) , but the objects proved to be too faint 
and/or without emission lines. The host of MACS0257 
cand2 is too faint for spectroscopy, with I > 27 mag. 
The luminosity of the SN candidate does not rule it out 
as a SN la. MACS0257 cand4 is offset from the core of 
a resolved host with an elliptical morphology. The host 
colors are consistent with the cluster red sequence in the 
available Subaru bands (/, z'). MACS1149 cand2 is off- 
set from the core of a host with an early-type morphol- 
ogy. The host colors are consistent with the red sequence 
in V — I, but are redder than cluster galaxies in R — I 
and R — z' . Spectroscopy of this galaxy failed due to a 
mistake in mask designs, as a result of which the galaxy 
spectrum feU off the CCD. MACS0647 cand4 is offset 
from the core of a resolved host with unclear morpho- 
logical type. Subaru colors indicate that it is ~ 1.5 mag 
bluer than the cluster red sequence. 

6.5. Classification Summary 

We conclude that five of the cluster events are likely 
SNe la, and one is likely a CC SN. To these we add, as 
possible cluster SNe la, the two hostless SNe, the two 
SNe with ambiguous hosts, and the four candidates for 
which we do not have measured redshifts. In terms of 
radial distribution, the cluster SNe are found at projected 
distances of up to 0.7 Mpc from the BCG (see Table S]). 

In order to reduce our uncertainties due to misclassi- 
fication, when calculating the SN la rate we impose a 
magnitude cut and remove from our sample candidates 
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fainter than 26 mag. This changes the SN count to five 
likely and six possible cluster SNe la. 

7. SN RATE CALCULATION 

With the sample of potential cluster SNe la con- 
structed above, we can now derive a SN rate for our 
cluster sample. The SN rate per unit stellar luminosity 
is calculated as follows: 

^la = , (1) 

j 

where N is the number of SNe, Aij is the visibility time 
(or "control time," the time during which a cluster SN la 
is above the detection limit of the jth image), i&andj" is 
the cluster luminosity within the search area of the jth 
image in the chosen photometric band, and the summa- 
tion is over all the survey images. The details of each el- 
ement in the calculation are given below. To account for 
the statistical nature of some of the quantities that enter 
the rate calculation, we conduct a Monte-Carlo simula- 
tion in which we measure the rate many times, each time 
assigning values to the measured quantities by drawing 
them from a distribution (see § 17. 3p . The final rate and 
its uncertainty is measured from a histogram of the re- 
sults. 

7.1. Visibility Time 

The visibility time depends on the detection efficiency, 
the peak luminosity, and the shape of the light curve at a 
given redshift and filter. In principle, it can also depend 
on the time interval between observations, as a SN is less 
likely to be detected via subtraction if it has a similar 
brightness at both epochs. In practice, since all of our 
comparisons are for epochs that are separated by at least 
one year, old SNe would have had enough time to decline 
below our detection limit, and thus the visibility time is 
not affected. We calculate the effective visibility time 
from 

/oo 
r][m{t)\dt, (2) 
-oo 

where m{t) is the SN la light curve in the image band- 
pass (either F814W or F775W) at the given redshift, and 
'rj[m{t)\ is the detection probability as a function of SN 
magnitude. We describe below each step in this calcula- 
tion. 

7.1.1. Detection Efficiency Estimate 

To determine the detection efficiency of our survey, we 
conducted efficiency simulations following the scheme de- 
tailed by Gal- Yam et al. (2002) and Sharon et al. (2007). 
To each field, we added some 200 fake SNe, in a range of 
magnitudes, and with a spatial distribution that follows 
the flux of the galaxies. The simulated images then un- 
derwent the same search procedure as the real data, and 
the number of SNe that were recovered in each magni- 
tude bin was noted. We find that, while the efficiency 
function strongly depends on the limiting magnitude of 
the image, its shape is also sensitive to other attributes 
of the field and of the observation. Although HST im- 
ages do not suffer from atmospheric distortion of the 
PSF, the cores of bright galaxies are not perfectly sub- 
tracted in difference images, resulting in residuals that 



can hide faint transients. This effect is amplified when 
two epochs are not obtained at the same position angle, 
due to position-dependent variations in the PSF. We ex- 
perimented with several PSF- matching techniques (e.g.. 
Gal- Yam et al. 2008) and determined that a simple sub- 
traction is sufficient for our purposes. We find that there 
is 100% efficiency in all epochs to detect SNe at magni- 
tudes brighter than 23, even in the cores of bright galax- 
ies. The efficiency differs from field to field, and even 
between epochs of the same field. Figure |6] shows the re- 
sults of our efficiency simulations. We parametrized our 
efficiency curves as a function of magnitude m with the 
function 

{^1 + e ir^\ , m < mo.5 
( 1 + e =2 j , m > mo.5, 

(3) 

where mo.5, s, and S2 are free parameters that are fit 
to the simulated efficiencies: mo.5 is the magnitude at 
which the efficiency drops to 0.5, and s and S2 determine 
the range of m over which 77 changes from 1 to 0.5 and 
from 0.5 to zero, respectively. 

In the final rate calculation, we apply a magnitude cut 
at 26 mag, by setting the efficiency above 26 mag to zero, 
and rejecting from the sample SNe that were detected 
above that value. 

7.1.2. Light Curves 

The light curves of SNe la exhibit an empirical inho- 
mogeneity, where luminous SNe tend to rise and decline 
more slowly than subluminous ones. This means that lu- 
minous SNe will be above the detection limit for longer 
than subluminous ones, resulting in an overall visibility 
time that is correlated with SN brightness. The cor- 
relation between peak magnitude and light-curve shape 
can be parametrized by a stretch relation (Phillips 1993; 
Perlmutter et al. 1999; see Leibundgut 2001 for a re- 
view), which exhibits an intrinsic root-mean square (rms) 
scatter of ~ 0.2 mag. In this paper, we use the form 
Ms = Ms=i—a{s — l), ts = ts^ixs, as described by Perl- 
mutter et al. (1999), where s is the stretch factor (s = 1 
means an unstretched light curve), Ms is the absolute 
i?-band magnitude, and a = 1.47, based on the empiri- 
cal results of Knop et al. (2003). The peak-magnitude 
distribution of SNe la also depends on the host-galaxy 
properties. By comparing light curves of SNe in different 
environments, Sullivan et al. (2006b) have shown that 
SNe la in elliptical galaxies tend to be dimmer (with a 
smaller stretch factor) than those in star-forming galax- 
ies (Hamuy et al. 2000; HoweU et al. 2001). 

For each cluster, we compile a set of light curves rele- 
vant to the cluster redshift and imaging band for several 
light-curve shapes and peak magnitudes. We start with 
a rest-frame, non-stretched, B-band template light curve 
from Nugent et al. (2002) and transform it to various 
stretched light curves, using the stretch relation above, 
and the measured distribution of elliptical host-galaxy 
stretch factors from Sullivan et al. (2006b). For consis- 
tency, we use the same non-stretched peak i?-band abso- 
lute magnitude as Sullivan et al. (2006a), Mb = —19.25 
mag (for Hq = 70 km Mpc~^). We assume an uncer- 
tainty of 0.15 mag in Mb, from the dispersion in peak 
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Fig. 6. — Detection efficiency curves for a point source in a differ- 
ence image, as a function of magnitude. The difference images are 
Epoch II vs. Epoch I (indicated by "cycle 14") and Epoch III vs. 
Epoch II (indicated by "cycle 15"). Magnitudes are measured in 
the difference image, and are Vega-based /-band magnitudes, ex- 
cept where indicated by "i." Cluster names are as marked. Circles 
mark the fraction of detected fake SNe, with error bars based on a 
binomial distribution. Solid lines are the best-fit efficiency curves 
with the indicated parameter values, parametrized as in Eq. [S] 

magnitudes of local SN la light curves after applying a 
stretch correction (Guy et al. 2005). This uncertainty is 
taken into account in our error budget (see ij 17. 3|) . 

A set of multi-epoch spectral templates from Nugent 
et al. (2002) are normalized to fit the S-band rest-frame 
stretched light curve, and rcdshiftcd according to the 
cluster redshift. We then combine the redshifted spec- 
tra with the HST bandpass in which the cluster was 
observed, to obtain a light curve for each combination 
of stretch, redshift, and bandpass in our survey. When 
calculating the visibility time for a particular image, we 
draw for each SN in the Monte-Carlo simulation a stretch 
factor (with its corresponding properly normalized light 
curve) from the Sullivan et al. (2006b) distribution of 
stretch factors. Table [T] lists the visibility time for each 
epoch, for the most likely stretch factor. We note that 
since epochs I and II are searched against each other, the 
values enumerated in Table [1] (i.e., visibility time, L, M, 
and the search area) are the same for epochs I and II and 
are only listed once. 

7.2. Cluster Stellar Luminosity 

A SN rate in a targeted galaxy population needs to 
be normalized by the stellar luminosity or mass of that 
population. These stellar properties, in turn, must be 
measured to some level of accuracy. Our rate measure- 
ment accuracy is limited by the small number of SNe that 



were discovered, with a la statistical error of order 30- 
40% (lower and upper Poisson errors for 10 events; see 
i) 17.3.11 below). An accuracy of 10% in the cluster stel- 
lar luminosities will therefore be more than satisfactory 
for our purposes. We generally follow the route detailed 
by Sharon et al. (2007), and perform "aperture pho- 
tometry" of the cluster light in several bands, as follows. 
For the MACS clusters and CLJ1226.9-H3332, we use the 
Subaru data centered on each cluster (see §[3]). For most 
of the MACS clusters, we were able to use calibrated V, 
R, and / images, with additional B- and z'-band imaging 
for some. 

Using the galaxy catalogs described in § [31 we measure 
the total flux in galaxies within a given aperture, and 
subtract from it the flux per unit area in galaxies in a 
"background" area at the outskirts of the images, multi- 
plied by the area of the aperture. This way, we statisti- 
cally subtract the contribution of "background" objects 
(those in front of and behind the cluster). The aper- 
ture in which we measure the cluster stellar luminosity 
is the exact search area (i.e., the area of overlap between 
epochs). The "background" is sampled in an annulus 
with area of at least 100 arcmin^, with an inner radius of 
at least 7' from the cluster center. These "background" 
annulus radii vary among the sample clusters, depend- 
ing on the angular size of the cluster and the available 
Subaru field of view, and take into account masking of 
bright foreground objects. The completeness of the Sub- 
aru data varies between clusters and depends on expo- 
sure time and observational conditions. For each cluster, 
we estimate the completeness magnitude as the turnover 
magnitude in its galaxy-magnitude histogram in the R 
band. We ignore in our calculation galaxies fainter than 
the completeness magnitude, and correct the total lumi- 
nosity accordingly as explained below. To avoid contam- 
ination from bright low-redshift objects, we also ignore 
galaxies that are brighter than the BCG. 

The luminosities of the four clusters for which we do 
not have Subaru data were measured in the same man- 
ner, but from the HST data. Since the ACS field of 
view is too small to select "background" annuli around 
the clusters, we assumed a universal background, drawn 
from the archival HST/ ACS imaging of The Great Ob- 
servatories Origins Deep Survey (GOODS; Giavalisco et 
al. 2004). 

For comparison with other SN rate measurements, we 
convert the observed net integrated galaxy-light fluxes 
in the available bands to rest-frame luminosities in sev- 
eral bands and form a cluster spectral energy distribution 
(SED). To each cluster SED, we fit redshifted template 
spectra of combinations of several types of galaxies, using 
synthetic photometry. For elliptical galaxies, we assume 
the Bruzual & Chariot (2003) synthesis of stellar pop- 
ulation with a single formation epoch at z = 3. Other 
galaxy types are represented by templates from Kinney 
et al. (1996). Typically, within each search area (i.e., 
the innermost ~ 1 Mpc of each cluster in our survey), 
the best fit was reached for a combination of 70-80% 
elliptical-galaxy fiux and 20-30% Sbc-galaxy flux. The 
fact that there is a non- negligible fraction of blue galaxies 
is consistent with recent measurements of the red fraction 
of galaxies in clusters in the redshift range of our sur- 
vey (Loh et al. 2008) due to the Butcher-Oemler (1978, 
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1984) effect. The luminosities in the desired rest-frame 
bands are measured via synthetic photometry of the best- 
fit template combination, scaled to fit the observed net 
cluster fiux. 

To account for undetected cluster galaxies, we correct 
the luminosity by multiplying it by the fraction of light 
that comes from the faint end of a Schechter (1976) lu- 
minosity function 



C 



S^mL)dL 



ilim("llim,z) 



(4) 



where ^{L)dL = $*(L/L*)"exp(-L/L*)d(L/L*). We 
adopt a = -1.00 ±0.06 and M* = -22.01 ±0.11 mag as 
the mean values for the (7-band luminosity function pa- 
rameters in clusters (Goto et al. 2002). We note that the 
depth of the Subaru images enables detection of galax- 
ies down to Ml = — 17 mag, ~ 5 mag fainter than M* , 
and thus the necessary correction is small, typically 1- 
5% for the lower- redshift clusters, and up to ~ 11% for 
the clusters at z > 0.8. 

Finally, to account for passive evolution of the domi- 
nant elliptical galaxy component over the redshift range 
of the cluster sample, the Bruzual & Chariot (2003) com- 
ponent is passively evolved backward or forward, as ap- 
propriate, to the mean visibility- weighted redshift of the 
sample, <z>= 0.6 (see Table [Ij. The Sbc component is 
not evolved in this calculation. This correction amounts 
to a 2% change in the total luminosity. 

7.3. Error Budget 

In this section we estimate the sources of uncertainty, 
both statistical and systematic. While the statistical er- 
rors can be propagated in a straightforward manner, the 
systematic errors affect the final result in a more compli- 
cated way, and may be correlated. To assess the overall 
systematic uncertainty, we calculate the SN rate by per- 
forming a Monte-Carlo simulation in which we measure 
the rate many times, each time assigning values to the 
measured quantities by drawing them from a distribu- 
tion centered on the best value, with a width according 
to the uncertainty of this value. Where applicable, we 
draw the values from a measured distribution and other- 
wise assume a Gaussian distribution. 

7.3.1. Statistical Errors 

The counting of SN explosions obeys Poisson statis- 
tics, from which we derive the statistical uncertainties. 
Contrary to the results presented for z w 0.2 clusters by 
Sharon et al. (2007), the number of secure cluster SNe 
in the present survey is not certain, due to uncertainty 
in SN classification and redshift, and is between 5 and 
11. We will consider this classification uncertainty as a 
systematic error below. If we adopt for the central num- 
ber of cluster SNe la the mean in this range, 8, the la 
Poisson errors are +49%, -35% (Gehrels 1986). 

7.3.2. Systematic Errors 

SN classification uncertainties. — A significant source of 
error in our rate derivation is the possible misclassifica- 
tion of SNe. Since our sample is not spectroscopically 
confirmed, some of the SN candidates may possibly be 
CC SNe at the cluster redshifts. CC SNe are prefer- 
entially located in star-forming regions, and are rarely 



found in the elliptical galaxies that dominate galaxy clus- 
ter environments. Although some cluster star formation 
may still be in progress, especially in high-redshift clus- 
ters, Saintonge et al. (2008) find that such activity tends 
to be outside the central 1 Mpc, while all of our cluster 
candidate SNe are within 0.7 Mpc in projection. Nev- 
ertheless, while many of the cluster-SN candidates in 
our sample occurred in elliptical galaxies, we note that 
some are associated with galaxies that have a late-type 
morphology. We will therefore conservatively allow for a 
maximal misclassification error. 

As described in detail in ^ we set a firm lower limit on 
the number of detected cluster SNe la, from the five can- 
didates in early- type cluster galaxies. Eight candidates 
are considered possible cluster SNe la, of which two are 
fainter than 26 mag and are rejected from our sample. 
This sets the upper limit due to classification uncertainty 
at 11. 

Luminosity error. — The derived cluster stellar luminosi- 
ties depend on several assumptions, such as the choice 
of area from which to draw the background, photometric 
errors, and the assumption of the galaxy templates with 
which the cluster SED was converted to rest-frame lumi- 
nosity. To account for these uncertainties, we perform 
a Monte-Carlo simulation, similar to those presented by 
Sharon et al. (2007). To assess the uncertainty intro- 
duced by our choice of background area, we explored dif- 
ferent choices of backgrounds for each cluster. We find 
that for a < 10% variation in the inner radius of the 
background annulus, the measured luminosity changes 
by ^ 5%. Skewing our measurement of the fraction of 
elliptical light in the cluster by 0.1 results in < 10% 
change in the measured luminosity. Assigning a 10% un- 
certainty to the luminosity of individual clusters in the 
Monte-Carlo simulation results in a 3.4% rms variation 
of the SN rate. 



Efficiency error. — As explained in i j7.1.1l we represent 
each efficiency curve by a continuous function that we 
fit to the results of efficiency simulations. To assess the 
errors induced by this process, and by the fact that the 
number of fake SNe from which the efficiency is derived is 
finite, we consider a range of efficiency-function parame- 
ters distributed normally around the best-fit value, with 
a 10% dispersion. In each realization in the Monte-Carlo 
simulation we draw a set of parameters from this distri- 
bution, resulting in an uncertainty due to the efficiency 
estimation of 1.6%. 

Visibility time error. — The visibility time depends on 
the shapes and peak magnitudes of the SN light curves at 
the cluster redshifts. These values are correlated through 
the stretch relation (see i i7.1.2|) . We explore a range of 
stretch parameters, following the distribution published 
by Sullivan et al. (2006b). In each realization in the 
Monte-Carlo simulation we draw a stretch factor from 
this distribution, and assign a light curve accordingly. 
Since the visibility time is not linear with a change of the 
stretch factor, the resulting distribution of SN rates is no 
longer centered on the value that is calculated from the 
most probable stretch factor. The inferred uncertainty 
due to the visibility time error is 2.7%. 

Accounting for all of the above systematic errors in 
the Monte-Carlo simulation, we get an overall system- 
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TABLE 5 
Cluster, SN Ia rate at <z>= 0.6 



Units" 


SN Rate 


Statistical 
Error 


Classification 
Error 


Systematic 
Error 


SNuM 
SNufl 
SNufl,o 


0.112 

0.35 

0.68 


+0.055 

-0.039 

+0.17 

-0.12 

+0.34 

-0.24 


±0.042 

±0.13 

±0.26 


±0.005 

±0.01 

±0.03 



^ SNus denotes SNe (100 yr IO^^Lb,©)" ■ SNum denotes SNe 
(100 yr IO^'^Mq)"^], and SNub.o denotes a rate normalized to 
cluster luminosity that is passively evolved z — 0. 

atic uncertainty of 5%. This uncertainty is, at present, 
negligible compared to the statistical errors, but will be- 
come relevant in future surveys that detect more than 
several hundreds of SNe. 



^ D.4 - 
0.35 - 

o 

"2 0.3 - 
S 0.25- 



0.2 - 



0.4 0.6 0.8 



7.4. Results 

The visibiUty-time-weighted mean redshift of our clus- 
ter sample is <z>— 0.6. For the adopted central value 
of iVia = 8, our measured SN rate per unit B-band 
stellar luminosity for this cluster sample is 0.35 SNub 
12 (statistical) ±0.13 (classification) ±0.01 (system- 
atic). 

The SN rate per unit luminosity at <z>— 0.6 cannot 
be compared easily to rates at low redshifts, even in the 
absence of any star formation, because of significant pas- 
sive luminosity evolution. For example, between z = 0.6 
and z = 0, a passive population fades by about a factor 
of 2 in blue bands (Bell et al. 2003; van der Wei et al. 
2005). To facilitate such a comparison, we also list in Ta- 
ble [5] a rate in units of SNub^q which is a cluster SN rate 
at <z>= 0.6 per unit stellar luminosity, but after pas- 
sively evolving forward, to z = 0, the elliptical Bruzual & 
Chariot (2003) component fit to each cluster's photome- 
try (the blue component, which contributes negligibly to 
the mass, is ignored in this calculation). 

To derive the rate normalized by stellar mass, we follow 
Mannucci et al. (2005), who converted iC-band galaxy 
luminosities to mass using the mass-to-light ratio derived 
by Bell et al. (2001) and observed galaxy B — K col- 
ors. We derive the rest-frame g — r colors of each clus- 
ter from the best-fit combination of template spectra. 
The stellar masses of the red and blue components of 
each cluster are then estimated from the color-dependent 
stellar mass-to-light ratio derived by Bell et al. (2003), 
logio(M0/Lg,o) = -0.499 + 1.519(5 - r) (see Mannucci 
et al. 2005 for a discussion of the validity of this ratio for 
our purpose). The total stellar mass in each search area 
is enumerated in Table [T] The SN rate is calculated from 
Eq. [1] in which we replace the stellar luminosity with the 
inferred total mass within the search area of each clus- 
ter. For A^ia = 8, the resulting SN la rate per unit stellar 
mass is 0.112 SNum to.039 (statistical) ±0.042 (classifi- 
cation) ±0.005 (systematic). The different rates derived 
in this work are listed in Table [5j 

8. DISCUSSION AND SUMMARY 

Our newly measured SN rate at z « 0.6 is the most 
accurate to date for clusters at such high redshifts. Com- 
bined with previous cluster-rate measurements, we are, 
for the first time, in a position of being able to examine 



Fig. 7. — Cluster SN rates from (in order of increasing redshift) 
Mannucci et al. (2008), Dilday et al. (2010), Sharon et al. (2007), 
Dilday et al. (2010), Gal- Yam et al. (2002), Graham et al. (2008), 
this work, and Gal- Yam et al. (2002). Triangles represent the 
SN rate derived from the upper and lower limits on the observed 
number of cluster SNe la, 11 and 5, respectively. Vertical error bars 
are Icr uncertainties, and horizontal error bars show the cluster 
samples' redshift ranges. 

the evolution of the SN rate in clusters over the greater 
part of cosmic history. Figure [7] shows the cluster SN 
rate as a function of redshift, at <z>= 0.02 (Mannucci 
et al. 2008), <z>= 0.08 and 0.22 (Dilday et al. 2010), 
<z>= 0.15 (Sharon et al. 2007), <z>= 0.25 and 0.90 
(Gal- Yam et al. 2002), <z>= 0.46 (Graham et al. 2008), 
and <z>= 0.6 (this work). For consistency, the rates 
from Gal- Yam et al. (2002) were converted from SNub 
to SNua/ using the ratio found by Sharon et al. (2007) 
for the lower redshift bin, and using the ratio found in 
this work for the higher redshift bin. The SN rate in 
the redshift bin studied in this work is consistent with 
the lower redshift rate measurements (to within uncer- 
tainties), and shows that there is, at most, only a slight 
increase of cluster SN rate with redshift. 

In a companion paper (Maoz et al. 2010b), we ana- 
lyze the results and discuss the implications of the ob- 
served rates on SN progenitor models, and on the role 
of SNe in the metal enrichment of the ICM. However, 
some conclusions emerge directly from the low SN rate 
observed out to z « 1. Independent of any model, the 
integral over the observed SN la rate per unit mass be- 
tween z = 1 and gives the total number of SNe per 
unit stellar mass. If multiplied further by the mean iron 
yield of a SN la, O.7M0 (e.g., Mazzah et al. 2007), one 
obtains M-pc/M^, the ratio of the iron mass produced in 
clusters over that cosmic period to the present-day stel- 
lar mass. Assuming, for example, a constant SN la rate 
of 0.1 SNum over the past 6 Gyr, as implied by Figure [71 
gives Mpc/M^, w 0.0004. In contrast, the observed ratio 
in present-day clusters, after subtracting the expected 
contribution from CC SNe, is [MFc/Af*]obs ~ 0.004 (see 
compilation and analysis of Maoz 2008), an order of mag- 
nitude higher. Thus, only a small fraction, < 10%, of the 
iron mass could have been produced by SNe-Ia between 
z = 1 and 0. This conclusion is strengthened if one as- 
sumes, instead of 0.7 Mq, the lower iron yields associated 
with subluminous SNe la. Such SNe tend to occur in the 
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early-type galaxies that dominate galaxy cluster cores. 
The factor ^ 2 increase in ICM iron abundance between 
these two redshifts, recently reported by Balestra et al. 
(2007) and Maughan et al. (2008), cannot be the ef- 
fect of new iron production during this time interval, and 
must instead be a redistribution effect. Alternatively, the 
abundance evolution may not be real (Ehlert & Ulmer 
2009), in which case the non-evolution is fully consistent 
with our SN results. 

Our measurements constrain the fraction of intergalac- 
tic SNe la (Gal- Yam et al. 2003). Assuming the two 
hostless SNe arc intergalactic cluster events, and that 
only the 6 most secure events in our sample are indeed 
SNe la in cluster galaxies, we get an upper limit on the 
relative fraction of intergalactic events of 2/8. This is 
very similar to the fraction measured by Gal- Yam et al. 
(2003) in lower redshift {z w 0.1) clusters (2/7). Since we 
cannot rule out that the two hostless sources we found 
are not cluster SNe la (i.e., that they reside in faint back- 
ground galaxies, or are not SNe at all), the lower limit on 
the intergalactic SN la fraction could be as low as zero. 
Assuming that 1 ± 1 of these events are cluster SNe la, 
and that our cluster SN la sample includes both the 6 
likely and 7 possible events we listed above, we get a 
likely fraction of 1/13 with an uncertainty of a factor of 
a few. This is several times below the value we measured 
at lower redshift (Gal- Yam et al. 2003), as would be ex- 
pected from models (Dubinski 1998) which predict that 
the intergalactic population of stars is mostly assembled 
between z = 1 and z = 0. 

To summarize, we have conducted a survey for SNe in 
the fields of 15 galaxy clusters, using new and archival 
HST/ ACS data. Each cluster was visited two to three 
times, for at least one orbit. In these data, we have dis- 
covered 37 candidate transient events, of which five are 
likely cluster SNe la, eight are possible cluster SNe la, 
and the rest are background or foreground events. We 
have determined cluster membership of the candidate 
host galaxies using follow-up spectroscopy from ground- 
based telescopes, and have measured stellar luminosities 
using Subaru and HST photometry. 

We find that the SN rate in clusters at 0.5 < z < 0.9, 
which we have measured here accurately for the first 
time, is consistent with the rates measured at lower red- 
shifts. Our main finding is thus that there is little or no 
evolution in cluster SN rates from the present time out to 
z w 0.9. The low and unevolving SN rate suggests that 
an increase in ICM iron abundance between redshift 
and 1, as reported based on X-ray observations, if real, 
is the result of the; rcidistribution of iron in clusters, and 
not due to the production of new iron by SNe during this 
period. Two of the candidate events are possible hostless 
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cluster SNe la, which we have discussed in the context of 
the few known examples of such intergalactic SNe, and 
their fraction in clusters. 

In a forthcoming paper, we will combine our current 
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